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AN ANALYTICAL STUDY OF THE EFFECTS OF JETS 
LOCATED MORE THAN ONE JET DIAMETER ABOVE 
A WING AT SUBSONIC SPEEDS 

By Lawrence E. Putnam 
Langley Research Center 


SUMMARY 


A procedure has been developed to calculate the effects of blowing 
two jets over a swept tapered wing at low subsonic speeds. The algorithm 
used is based on a vortex lattice representation of the wing lifting surface 
and a line sink-source distribution to simulate the effects of the jet 
exhaust on the wing lift and drag. The method is limited to those cases 
where the jet exhaust does not intersect or wash the wing. The predictions 
of this relatively simple procedure are in remarkably good agreement with 
experimentally measured interference lift and interference induced drag. 

The results of the analytical study have indicated that substantial 
Increases in wing lift and reductions in wing induced drag can occur 
when jets are blown over a wing. The magnitude of these effects increase 
with increasing ratio of jet exit velocity to free-stream velocity. The 
results also indicate that the interference effects are the most favorable 
when the nozzle exits are located ahead of or at the wing leading edge and 
inboard near the wing center line. The present method does not correctly 
predict the effects of jet vertical location at heights less than 
approximately 1.5 jet exit diameters; however, as the jet vertical location 
is reduced to this height the favorable effects of blowing jets over a 
wing increase. 
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INTRODUCTION 


Recently, there has been renewed interest in locating the propulsion 
systems of future jet transport airplanes forward and above the wing. This 
unconventional engine location can take two forms. In one form, the jet is 
exhausted tangential to the wing upper surface to take advantage of Coanda 
turning of the jet exhaust as it passes over the wing to achieve high lift 
(refs. 1 through 4). In the second form the jet exhaust is raised so that 
the exhaust flow does not attach to the wing to avoid the cruise drag 
penalties associated with the first form. In both cases the jet noise 
propagated downward during takeoff and landing may be reduced by the shield- 
ing effects of the wing (refs . 5 and 6) . In addition to noise shielding 
benefits, locating the propulsion system forward may also alleviate mass 
balance problems associated with advanced supersonic transport airplane 
concepts and may Improve the wing flutter characteristics of such airplanes. 

To date, however, except for references 7 through 9, there is very little 
information available for use in assessing the effects of over wing jet 
blowing on subsonic or supersonic cruise performance and on climb performance 
of future jet transport airplanes. 

The purpose of the present study was, therefore, to investigate analy- 
tically the effects of blowing jets over wings on subsonic cruise and climb 
performance. An elementary analytical procedure based on the vortex lattice 
theory for the wing and the theory of reference 10 for a jet exhausting into 
a subsonic stream has been developed to calculate the effects of blowing one 
or two jets over a swept tapered wing. As a result of the basic assumptions, 
the method is limited to the second case where the jet exhaust does not inter- 
sect or wash the wing. As long as the jet is sufficiently high such that it 
does not violate the above condition the effects of jet to free-stream 
velocity ratio (i.e. jet thrust coefficient), location of the jet exhaust 

relative to the wing, nozzle exit diameter, wing leading edge sweep, wing 
aspect ratio, and wing taper ratio on the lift and drag induced by blowing 

jets over a wing at subsonic speeds can be analytically investigated. 
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SYMBOLS 


Ae 


exit area per nozzle 


AR 


wing aspect ratio, 


ref 


wing span 


axial force coefficient 


T>»i 


induced drag coefficient 


AC 


D,i 


(C D ± ) jet on - (C^) jet off 


drag coefficient at zero lift 


lift coefficient 


AC, 


(C L ) jet on - (C L ) jet off 




normal force coefficient 
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c 


wing chord, or chord of elemental vortex lattice panel 


c 

1 


section lift coefficient 


& nozzle exit diameter 

e 


K, 


slope of line source or line sink distribution 


length of chordvlse vortex-line segment on wing 


l core 


length of core of jet exhaust plume 


M 


Mach number of jet at nozzle exit 




free stream Mach number 


N 


number of points or number of elemental vortex lattice 
panels 


‘t ,e 


total pressure of jet at nozzle exit 


free stream static pressure 


Q 


strength of line source or sink distribution 



dynamic pressure of Jet at nozzle exit 




free stream dynamic pressure 


ref 


reference area (wing planform area) 


semi-width of horseshoe vortex 


U g velocity of Jet at nozzle exit 


U 


free stream, velocity 



circulation induced perturbation-velocity 
components in positive X-, Y-, and Z- directions 


U jet* v jet’ w jet 


Jet induced perturbation-velocity components 
in positive X-, Y-, and Z- directions 


V 

e 


V 

r 


effective velocity ratio 



radial velocity in Jet axis system induced by 
line sink-source distribution 
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orthogonal right-handed primary Cartesian 
coordinate system with origin at wing apex 
(see figure l). Positive X-direction is 
forward, positive Y-direction is toward 
right wing tip, and positive Z-direction is 
downward. X-Y plane parallel to wing chord 
plane 


jet cylinderical coordinate system with origin 
at nozzle exit (figure 3); positive X-direction 
is in direction of jet exhaust 


cylinderical coordinates associated with jet 


longitudinal location in primary Cartesian 
coordinate system of midspan of quarter-chord 
of elemental panel 


longitudinal location in primary Cartesian 
coordinate system of midspan of three-quarter- 
chord of elemental panel 


primary Cartesian coordinates of nozzle exit 


primary Cartesian X-coordinate of wing leading 
edge at 


primary Cartesian Y-coordinate of control 
point on elemental panel 



angle of attack, deg 


circulation strength 


ratio of specific heats of jet at nozzle exit 


ratio of specific heats of free stream 


dummy variable of integration 


fraction of wing semispan from wing root chord in 
Y-directi.on 


location of jet axis in fraction of wing semispan. 



b 


camber angle of wing at control point, deg 


sweep angle , deg 


wing quarter chord sweep angle, deg 


wing leading edge sweep angle 



X 


taper ratio, 



p density of jet at nozzle exit 


o free stream density 

00 


Subscripts 


av 


c 


r 


s 


t 

total 


average 

chordwise 

root 

spanwise 

tip 

total 


METHOD OF ANALYSIS 

The algorithm used in the present study to calculate the effects of 
blowing exhaust jets over a swept tapered wing is based on a vortex lattice 
representation of the wing lifting surface and a line source-sink -distri- 
bution to represent the effects of the exhaust jets. The method, which has 
been programmed for a digital computer, assumes that the flow external to 
the jet exhaust is steady, irrotational, inviscid, and incompressible. It 
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was also assumed "that the jet is not deflected hy the free stream* the jets 
does not intersect or wash the wing, and the jet cross sectional shape is 
not distorted by the wing flow field or by any crossflow components of the 
free stream velocity. A typical wing planform and engine arrangement 
for which the effects of blowing jets over the wing can be calculated 
with the present method is shown on figure 1. 

The vortex lattice method used to represent the wing lifting 
surface is based on references 11 and 12. Each half of the wing is 
subdivided in both the spanwise and chordwise directions into 50 or 
less elemental areas. In the present computer program the number of 
spanwise and number of chordwise panels is arbitrary; however, unless 
otherwise specified, all calculations for the present study were made 
with 5 equal spaced chordwise and 10 equal spaced spanwise divisions 
on each half of the wing. Each elemental area is represented by a 
horseshoe vortex with the bound portion lying along the local 
quarter-chord-line of the element. (See figure 2.) The trailing 
vortices lie streamwise along the inboard and outboard edges of each 
panel on the wing chord plane. Aft of the wing trailing edge, the 
trailing vortices continue in the chordwise direction, that is, it 
is assumed that the trailing vortices are not turned by the wing 
downwash or by the free stream velocity. Note that this is not the 
conventional method of locating the trailing vortices parallel to 
the free stream velocity as used in reference 12, but is similar to 
the formulation of reference 11. The boundary condition that the 
flow be tangential to the wing surface is satisfied for each element 
at a point on the lateral mid-point of the local three-quarter-chord 
line of the element. 

The induced effects of the jet exhaust were simulated with a 
distribution of line sinks and sources located on the longitudinal 
axis of the jet. The axis of the jet was divided into a number of 
segments over which the sink or source strength was assumed to be 
linear. This continuous line sink-source distribution is equivalent 
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"to a. series of triangular elemental distributions as illustrated in 
figure 3. The radial velocity induced by the jet at a point in the 
surrounding flow field is therefore given by the following equation: 


V r <x,r) 



(x 2 - O d£ 

[U-O 2 * r 2 j 


3/2 




The strength of each elemental line source or sink was adjusted to 
give the inflow velocity on the boundary of the jet predicted by the 
method of reference 10. In the present analysis only those 
interference effects resulting from the jet radial perturbation 
velocities are considered: any perturbation velocities induced 

parallel to the jet axis are neglected. In effect, therefore, the 
interference effects of the Jet on the wing are attributed to an 
upwash distribution along the wing lifting surface. 



To determine the circulation strength of each elemental horseshoe 
vortex on the wing, the following equation for the circulation 
induced downwash velocity must be satisfied at each control point on the 
wing: 


V i = sin (a + 0) - ".let - “jet 


v . 


u. 


U 


cos 0 


U 


U 


tan 0 


A sketch illustrating this boundary condition at each control point is 

shown on figure U. It is apparent from this sketch how the induced effects 

of the jet modify the boundary condition equation. Note that in the 

present formulation the u perturbation velocity only arises when the 

j et 

axis of the jet is not parallel to the wing chord plane. 

After the circulation strengths for the horseshoe vortices are 
determined, the forces acting on the wing can be determined. From 
reference 11, the forces acting on each spanwise segment of vortex filament 
in coefficient form are: 


"A,s 


2r 

c 

~U~ 


w. 

l 

U 


u 



- 

1 

2s 


S . 
ref J 


•'N.s 


2r 

c 

u~~ 


fv. v, . \ u, u. , 

^ + r Jet ) t A _ J. _ + 


U 


u 


/ 


u 


u 


cos a 


2s 


ref 


For each chordwise segment of vortex filament the force coefficients are: 


C A =0 
A,c 


11 



N,c 


2T 

c 

U 


V. V . 

_i + J.fA 


U 


U 


ref 


Where T is the net circulation strength resulting from the 
c 

individual circulations of each trailing horseshoe vortex leg forming 
that segment. The total normal force and axial force acting on the 
wing are obtained by summing the normal and axial forces acting on all 
segments. The lift and induced drag coefficients are then obtained from: 


C T * C„ . * , cos a - C. , sin a 
Xj A j to lqJ. 


C = C cos a + CL , . , sin a 

D,i A, total N , total 


This solution for the induced drag is equivalent to the near field 
solution of reference 12. The integration of the wing span load 
distribution (the far field method) to obtain the induced drag is not 
applicable to the present problem. Because of the induced upwash field 
generated by the jet exhaust the assumption inherent in the far field 
method are violated. (See refs. 13 and 1*+.) 

DISCUSSION 

Comparison of Present Method with Other Theoretical Methods and 

with Experiment 

To verify the vortex lattice algorithm used to represent the wing 
in the present method, the computed lift and drag characteristics of an 
aspect ratio 8 straight wing without a jet have been compared to the 
predictions of reference 12. (See fig. 5*) The present computer 
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algorithm, gives essentially the same results as reference 12 for each 
of the cases considered: a flat wing, a wing with Clark Y camber 

distribution, and a twisted wing with no camber. There are some small 
discrepancies in the results of the two methods at the higher angles of 
attack. This small discrepancy is associated with the nonlinearization 
of the boundary conditions in the present method. Note that the near 
field induced drag results of reference 12 are compared with the 
present method. 

Experimental investigations of the effects of blowing jets over a 
wing are very scarce; in particular for the case where the jet is 
sufficiently high such that it does not wash the wing. One source of 
such data is reference 15 • The comparison of the present method with 
these experimental results shown on figure 6 indicates very good 
agreement at vertical locations of the jet exhaust greater than 1.5 
nozzle exit diameters above the wing. In general, the effects of jet to 
free-stream velocity ratio, longitudinal location of the jet, and vertical 
distance of the jet above the wing on the interference lift are well 
predicted. As the jet approaches the wing, however, the experimental 
results show a significant increase in interference lift over the 
predicted values. This discrepancy results when the jet washes the wing 
and Coanda turning of the jet with its associated increase in 
interference lift occurs. This effect of Coanda turning of the jet 
exhaust is of course neglected in the present method. 

A comparison of the present method with the experimental data of 
reference 9 for two Jets blowing over a 50° swept leading edge wing 
with an aspect ratio of 3 and a taper ratio of 0.3 is shown on figure 7. 
These experimental data were obtained at nozzle exit total pressure to 
free-stream static pressure ratios from approximately 1 to 6 and at free- 
stream Mach numbers of 0.U, 0.6, and 0.7. The theory of reference 10 was 
derived for an incompressible jet exhausting into an incompressible 
external stream with no density difference between the two streams. The 
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variable density case can be related to the incompressible constant 
density solution of reference 10 with an effective velocity ratio (ref. 1 6) 
as defined in the following equation 



For an air Jet with = Y^, - 



As shown on figure 7 (a) and 7 (b) the effective velocity ratio 
satisfactorily correlates the interference lift and the interference 
induced drag resulting from blowing two jets over the 50° swept leading 
edge wing of reference 9 at Mach numbers from 0.4 to 0.7. The predictions 
of the present method are in good agreement with the experimental results. 
These results indicate that the present method can be extended to the 
case where the Jet and free stream are compressible and at different 
densities by use of the effective velocity ratio concept. An illustration 
of the application of the present method using the effective velocity 
ratio concept to the prediction of the effects on lift and induced drag 
due to blowing jets over a wing is shown on figure 7 (c) for the 50° swept 
wing of reference 9. The following equations were used to obtain the 
predicted lift and induced drag when the jets were blowing over the wing: 
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C = (C ) „ + AC t .. , , at a given a 

L, predicted ' L,jet off measured L, predicted 


(r ) * lie ) ) + ( AC , ) , at a given C 

^D,i 'predicted ^ D,i jet off 'measured D,i predicted L 


The predicted lift and drag characteristics using the present method to 
calculate the effects of the jets with the above equations are in good 
agreement with the experimental results. Of particular interest is the 
experimentally indicated improvement in induced drag (not necessarily 
in total drag) due to blowing two j.ets over the wing and the similar 
predictions of the present method. 

Effect of Vortex Lattice Arrangement 

Shown on figure 8 is the effect of the number of elemental chordwise 
panels on the predicted lift and induced drag of an aspect ratio 8, 
taper ratio 0.3 wing having the quarter chord swept 0° with and without 
two jets blowing over the wing. It is apparent that there is essentially 
no effect of varying the number of elemental panels on the predicted 
lift and induced drag of the wing or on the increments in lift and drag 
due to Jet blowing. Varying the number of spanvise elemental panels 
does result in a change in the predicted lift and induced drag for the 
wing with and without jet blowing. (See fig. 9.) It appears, however, 
that the interference effect due to jet bloving is not a strong function 
of the number of elemental vortex lattice panels. Therefore, regardless 
of the vortex lattice arrangement used, the predicted effects of jet 
blowing are essentially uneffected. Therefore 5 chordwise and 10 spanvise 
elemental panels have been used for all calculations of the present paper. 
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Effect of Velocity Batio 


The predicted effects of effective velocity ratio on the 
interference lift and interference induced drag due to jet blowing 
are shown on figures 6, 7» s^d 10 for sb aspect ratio 2 straight wing, 
an aspect ratio 3, taper ratio 0.3 wing with 50° swept leading edge, 
and an aspect ratio 8, taper ratio 0.3 wing with 30° swept quarter 
chord, respectively. The present method predicts an increase in inter- 
ference lift coefficient with increasing Jet blowing (that is, 
increasing 1/V e ). This favorable interference effect is essentially 
uneffected by angle of attack. The present method also predicts a 
decrease in induced drag with increasing ratio of jet velocity to 
free-stream velocity l/V^. This favorable increment in interference 

induced drag increases with lift coefficient. On figure 10 (b) are 

A 

shown typical calculations of the effects of effective velocity ratio 
on wing span load distribution and on section lift coefficient. 
Increased jet blowing results in an increased distortion of the span 
load distribution. Note that this increased distortion does not imply 
increased induced drag since integration of the span load distribution 
to obtain induced drag is only valid for planar wings without external 
interference effects. As would be expected the span load does increase 
in the vicinity of the jet. Increasing jet blowing also increases 
the section lift coefficient over the entire wing span. 

Effect of Jet Location 

The predicted effects of jet location on the interference 
lift and interference induced drag due to two Jets blowing over 
various wing configurations are shown on figures 6 and 11 through l4. 
Increasing the vertical distance of the jets above the wing (figs. 6 
and 11 ) causes a decrease in the Interference lift and a reduction 
in the favorable interference induced drag. Moving the nozzle exit 
ahead of the wing leading edge causes a small increase in 
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interference lift and a small increase in the favorable induced drag 
increment. (See fig. 12.) Moving the nozzle exit aft of the ving 
'‘’leading edge causes a decrease in the interference lift and a 
substantial reduction in the favorable induced drag increment. The 
calculations shown on figures 13 and lh indicate that the closer the 
jets are located to the wing center line the more favorable the 
interference effects due to blowing two jets over a wing. This is 
indeed a fortuitous result since reference 4 indicates that for a 
two engine configuration the engine should be located inboard to 
minimize the engine-out effect on rolling moment associated with upper 
surface blowing configurations during takeoff and landing. 

CONCLUDING REMARKS 

A procedure has been developed to calculate the effects of 
blowing two jets over a swept tapered wing at low subsonic speeds. 

The algorithm used is based on a vortex lattice representation 
6 of the wing lifting surface and a line sink-source distribution 
to simulate the effects of the Jet exhaust on the wing lift and 
drag. The method is limited to those cases where the jet exhaust 
does not intersect or wash the wing. The predictions of this 
relatively simple procedure are in remarkably good agreement 
with experimentally measured interference lift and interference 
induced drag. 

The results of the analytical study have indicated that 
substantial increases in wing lift and reductions in wing induced 
drag can occur when Jets are blown over a wing. The magnitude 
of these effects increase with increasing ratio of jet exit 
velocity to free-stream velocity. The results also indicate 
that the interference effects are the most favorable when the 
nozzle exits are located ahead of or at the wing leading edge and 
inboard near the wing ceter line. The present method does not 
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correctly predict the effects of jet vertical location at heights 
less than approximately 1*5 jrt exit diameters; however » as the 
jet vertical location is reduced to this height the favorable 
effects of blowing jets over a wing increase. 
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Figure 2,- Sketch of e typical elemental vortex -lattice panel on wing, 







(a) Flat wing. 


Figure $♦- Comparison of present method with that of reference 12 -for an AR * 8, 

X =* 1.0 straight wing (jet off) * (Reference 12 near field solution for 
induced drag shown,) 





(b) Wing with Clerk Y camber distrib 


Figure 5.- Continued 





(e) Wing with twist distribution 


Figure 5.- Concluded. 






Figure 6. - Comparision of present method with experimental data of reference 15 for blowing one jet 
over an AR ■ 2, X - 1 straight wing. " 0. 192, (x fi - x^)/d e - 1.0, - 0.0, 

and a ■ 0°. 



AC|_ 


Figure 7.- 



2 3 4 5 


'/ v e 

(a) Interference lift coefficient* 

Comparison of present method with experimental data of 
reference 9 for two jets blowing over a wing with AR - 3, 

X ** 0.3, and = 50*. ** 0.0078, (* e - = 2.59, 

n = 0.46 , z /d = -1*5. 
e e e 
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(c) Effect of jet blowing on lift and drag coefficient at 


Figure 7. - Concluded. 









Figure 8. - Effect of number of chordwise elemental vortex lattice panels when N $ ■ 10 on the calculated lift and 
induced drag of a wing with AR ■ 8, A • 0.3, and ■ 0° and with two jets having A g / S re f. ” 0. 02; 
located at (x - x, )/ d - 0. 0, n ■ 0. 3, and ij d g ■ -1. 5. 
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Figure 9. * Effect of number of spanwise elemental vortex lattice panels when * 5 on the calculated lift and 

and induced drag of a wing with AR * 8, X * 0.3, and A^ 4 * 0® and with two jets having A g /S ref • 0.025 
located at (x g - x Ie >/d e * 0.0, • 0.3, and zjd e ■ -1.5. 
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(b) Effect, on span load distribution and 
section lift coefficient at a * 2°« 


Figure 10.- Continued. 












Figure 11,'- Calculated effect of vertical location of two jete on aerodynamic characteristics 
of an AR •* 8, ^ - 1,0 straight wing; A / s ef ~ 0,0125, V = 0.20, n e = 0.25, 







z e/ d e C L C D,i 


Jet off 0J647 0.0011 



(b) Effect on span load distribution and 
section lift coefficient at a ^2°. 


Figure 11.- Continued. 










(a) Effect on and 


Figure 12.- Calculated effect of longitudinal location of two jets on aerodynamic 

characteristics of anAR* 8, X** 1.0 straight wing. -^ e ^ re f * 

V 0 - 0,20, ti c = 0,25* z e / d e “ ^ 


*fo 


■1.401 



( x e — *Ze)/ d e C L c D,i 


Jet off 0.1647 0.001 l 

2.803 .2398 -0,0019 



(b) Effect on spaa load distribution and 
section 1 lift coefficient at a ■ 2°. 


Figure 12.- Continued. 
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(c) Effect on AC and AC. . . 


n 


Figure 12.- Concluded 



(a) Effect on and 

Figure 13,- Calculated effect of lateral location of two jets on aerodynamic 
characteristics of an AR a 8, \ *» 1.0 straight wing, A e ^ S re f ° 

V « 0.20, (x - x 0 )/d = 0.0, z /d - -1.401. Wl ' 
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Figure 13.- Continued. 



(c) Effect/ on ACL and AC. 


Figure 13. - Concluded. 
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